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Recently topological superconducting states has attracted a lot of interest. In this work, we consider a topo¬ 
logical superconductor with Z 2 topological mirror order (Tl and s±-wave superconducting pairing symmetry, 
within a two-orbital model originally designed for iron-based superconductivity m We predict the existence 
of gapless edge states. We also study the local electronic structure around an adsorbed interstitial magnetic 
impurity in the system, and find the existence of low-energy in-gap bound states even with a weak spin polar¬ 
ization on the impurity. We also discuss the relevance of our results to the recent STM experiment on Fe(Te,Se) 
compound with adsorbed Fe impurity O, for which our density functional calculations show the Fe impurity is 
spin polarized. 

PACS numbers: 73.20.r, 71.70.Ej, 71.10.Pm, 74.20.Pq 


Topological superconductor (TSC) has been attracted lots 
of interest on its potential application to the fault-tolerance 
quantum computation ID and superconducting spintron- 
ics El. The TS can be designed through the interface of a 
topological insulator with the fully gapped superconductor O 
or by engineering through the interface of ferromagnetic¬ 
superconducting nanostructure under several special condi¬ 
tions CHni. Most of these efforts were aimed to look for 
the edge states and in-gap bound states in a fully gapped su¬ 
perconductor. There are two types of edge states for a two- 
dimensional (2D) TSC. The first type is a chiral edge state 
in a time-reversal-symmetry-broken superconducting pairing 
state like [px + Wy\yy CHI [131 • The second type is a heli¬ 
cal edge state in a time-reversal-symmetry-invariant pairing 
state like [px ± Wy\yi) EdU. Recently, the crystalline sym¬ 
metry protected topological phase with mirror Chern number 
has been studied in a triplet pairing state of Sr 2 Ru 04 (TEh , for 
which chiral edge states are obtained M- A realization of 
helical edge state can be established in a junction formed by 
an s-wave superconductor and a strong topological insulator 
(TI), mimicking a time-reversal-symmetry-preserved p-wave 
superconducting pairing state 0. 

In this work, we propose a feasible scenario to realize a he¬ 
lical topological mirror superconductor (HTMS) based on the 
coexistence of the interaction-driven Z 2 topological mirror or¬ 
der and 5±-wave pairing symmetry. Our concept is based on 
the nature of a topological metal, where the superconductivity 
can be generated through a finite Fermi surface with electron- 
(hole-) Fermi pockets. The Z 2 mirror topological order cre¬ 
ates robust edge states even in the presence of a fully gapped 
5±-wave pairing symmetry. We further analyze the local elec¬ 
tronic structure of an adsorbed interstitial magnetic impurity 
(IMI) in the 2D bulk of the HTMS, and identify the different 
consequence from the topological and non-topological super¬ 
conductor. It is important to note that, since the supercon¬ 
ductivity is intrinsically inherited from the normal-state band 
structure, our proposal avoids the need to make a TI-SC het¬ 
erostructure for helical edge states. To make our model into a 
real-material context, we also carry out the DFT calculations 


for an adsorbed interstitial-Fe impurity in recently discovered 
11-family of iron-based superconductors FeSe, and find the 
adsorbed Fe-atom to be spin polarized. Therefore, our re¬ 
sults have a direct relevance to the recent STM experiment on 
Fe(Te,Se) SCs, where a robust zero-energy bound state (ZBS) 
has been found around an adsorbed Fe atom Cl. 

We consider the following effective Hamiltonian to de¬ 
scribe the coexistence of the Z 2 topological mirror order with 
the 5 ±-pairing symmetry, 

= T[h.6,n] + v[XAOH] + r[A]. ( 1 ) 

Here T is the kinetic energy with optimized hopping terms 
(ti-e) and chemical potential (p) to describe the low en¬ 
ergy physics that originally designed for the Fe-based com¬ 
pound O [TtI, V is responsible for the emergent Z 2 topo¬ 
logical mirror order Q with tunable anomalous orbital 
Hall order (Xaoh) for a (i-wave form factor As(k) = 
XAOH[cos{kx) - cos{ky)]a,a (whcrc Of G dxz 5 dyz, 
s Gti), and V is the s±-wave pairing order (A) with the form 
factor A(k) = 4A[cos(/ca^) x cos(/cy)]>t^ 4 , (in 1-Fe BZ). See 
definitions of T, V and V in the Method section. Through¬ 
out the work, the chemical potential is adjusted to ensure the 
system is half filled. In the following section, we numerically 
solve the band structure of Eq. O to reveal the interplay of 
topological (Xaoh) and superconducting (A) orders. 


I. EDGE STATES WITH FULLY GAPPED S± -WAVE 
PAIRING 

The TSC property can be revealed by studying the nature 
of edge states. We first calculate the energy dispersion un¬ 
der several conditions. Figure shows the evolution of the 
non-superconducting (NSC, A = 0) band structure with dif¬ 
ferent values of Xaoh in Brillouin zone (BZ) corresponding 
to the 2-site per unit cell. In Fig. [^, the band crossing be¬ 
low the Fermi energy occurs along the M-T direction when 
there is no existence of the anomalous orbital order (AOH), 
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FIG. 1: The band structures shows the gapless edge states of the coexistence of topological and superconducting orders, a, The folded 
band structure of a 2-site per unit cell BZ. b-e, The band structure of a strip of width of 20 lattice constants with open and close boundary 
condition under non-superconducting and superconducting phase. We take the periodicity along the ^-direction of the strip, where 2000 ky- 
points are taken. The inset of the middle panel of e is the enlargement of the crossing point of the Fermi surface, where we can find that there 
are totally four degnerated points around ky — (the green inset) and the others around ky — ^ and zb tt are not degenerated. 
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FIG. 2: Cartoon picture for an interstitial impurity in a 2D lat¬ 
tice. a. The top-view of the parent lattice structure and the orientation 
of dxz and dyz orbitals, b. The position of an interstitial atom located 
at the center of a plaquette of the 2D square lattice together with the 
local hopping process among the dxzidyz orbitals, where tc 7 , 7 r stand 
for orbital lobes aligned parallel (perpendicular) to the hopping di¬ 
rection between two nearest-neighboring sites of the square lattice). 
We set ter = atn = ^ for the parameterization. 


^AOH = 0, and electron and hole pockets are formed at the 
Fermi energy. In the presence of a small value of the AOH 
term (Xaoh = 0.15), the band crossing is lifted and the four 
bands are formed into two disentangled groups. However, the 
Fermi pockets remain. Furthermore, when the topological or¬ 


der is strong enough (Fig.[^ with Xaoh = 0.6), an indirect 
band gap is open directly at the Fermi energy and the bulk 
system becomes insulating, from which no superconducting 
pairing can emerge. Now, we construct a strip geometry to 
study the electronic states at the edge with periodic (Fig.[^,d) 
and open (Fig.[T]:,e) boundary conditions. Since there still ex¬ 
ists the translational invariance along the strip direction (that 
is, ^-direction) at the original lattice constant, the momen¬ 
tum along this direction is a good quantum number, which 
allows to solve the problem for each given wave vector ky in¬ 
dependently. We can see that the edge bands show up by com¬ 
paring Fig. [^,d and c,e. In the metallic/topological-metallic 
states(AAOi^ = 0, 0.15) without the superconducting order, 
the edge bands have curved features across the Fermi surface 
as highlighted in red in Fig. [T]:. For Xaoh = 0 the non¬ 
superconducting (NSC) edge states are four-fold degenerate 
(including two spin orientations and two edges) at a given ky 
near the region of ky = When Xaoh 0, the degen¬ 
eracy of edge bands near ky = ±f is reduced to two fold 
(that is, only two-fold spin degeneracy) at a given ky. How¬ 
ever, these bands are still crossing the Fermi energy. We now 
turn on the s±-wave pairing (A = 0.02) (see Fig.[^,e). In 
the absence of the topological order, the edge states have no 
zero-energy modes, which means the quasiparticle states are 
always gapped. Only in the presence of a sufficiently large 
topological order, can the edge states appear to cross the Fermi 
energy (see Fig.[J with Xaoh = 0.15, 0.6). We note that 
for Xaoh = 0.15, the normal state is metallic and the super¬ 
conducting pairing is physically feasible. Therefore, we will 
focus on the case of Xaoh = 0.15 in the next section. For 
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FIG. 3: LDOS and band dispersion from the model Hamiltonian, a, The calculated LDOS for various values of spin exchange interaction, 
J, on the impurity site (blue lines) and a site deep into the bulk (red lines) as a reference, b, The highly folded band structure for a supercell (of 
size N X N) lattice with a single interstitial impurity sitting in the center of each supercell, c, The LDOS on the interstitial impurity site (red- 
solid lines) and a site deep in the bulk (black-solid lines) for Xaoh = 0.0 (top panel) and Xaoh = 0.15 (bottom panel) at a fixed J = 0.8. 
d, The spatial dependence of LDOS at zero energy in the entire supercell of the size 50 x 50. The inset is the cut from site indeces (-24,0) 
to (25,0) along a bond direction through the one side of the plaquette, in which the interstitial impurity sits (with site coordinate (0.5,0.5)). 
Unless specified otherwise, the order parameter values are Xaoh = 0.15 and A = 0.02. 



bound state follows closely follow the symmetry of that of 
Majorana fermions. Inspired by the observation of a robust 
zero-bias conductance peak in the STM experiments around 
an interstitial-Fe-impurity on the surface of the Fe(Te,Se) 
compound a, we construct a microscopic model for the lo¬ 
cal adsorbed impurity within our model. Figure shows the 
2D lattice structure for a pristine system, where -orbital 
orientation is considered; while Fig. shows an adsorbed 
Fe atom and its effective local hopping with the neighboring 
Fe atoms in the parent compound. We write down the total 
Hamiltonian with the local effect of absorbed site, 


FIG. 4: Band structure from the model Hamiltonian, a, Xaoh = 

0 and b, Xaoh = 0.15. The high-symmetry momentum points of 
r, X" and M" are defined according to the Brillouin zone of a 1- 
site per unit cell. The superconducting order parameter A = 0.02 is 
taken. 

^AOH = 0.6 (see third row of Fig. [^), the normal state is 
already insulating and the superconducting pairing is not pos¬ 
sible. However, we include this case to show the structure 
of edge states in response to a “superconducting” particle- 
particle pair. 


II. LOCAL ELECTRONIC STRUCTURE AROUND AN 
ADSORBED IMPURITY 

The TSC property can also manifest in the local electronic 
structure around a single impurity. The use of a single impu¬ 
rity is a powerful method to probe the superconducting pairing 
symmetry Gll. More recently, a universal impurity-induced 
mid-gap state has also been predicted as a clear signature of 
topological superfluids |[T9l [20]|, where the near-zero-energy 


^Tot ^ yTSC ^ Jimp [ ^ ^ gimp j j ] ^ (2) 

in which is the local hopping term from the impurity 

site to its nearest neighboring Fe atoms in the parent sys¬ 
tem with tunable hopping parameter ^ (see Fig.]^), 
describes the local spin-polarized energy levels on the impu¬ 
rity site with an exchange parameter J. See more details for 
j-imp gimp Methods Section. We will focus on 

the local density of states (LDOS) for fixed topological order 
^aoh{= 0.15 and 0.0) and SC pairing strength A(= 0.02) 
with varying ^ and J. 

In the regime of a weak hybridization, one can anticipate 
that the resulting peaks in the LDOS on the adsorbed impu¬ 
rity site merely reflects the energy levels of an isolated atom. 
In our work, we are more interested in the regime where the 
coupling between the adsorbed impurity and its neighbors is 
relatively strong. Figure shows the LDOS calculated on 
the adsorbed impurity for various values of ^ and J but with 
fixed A = 0.02 and Xaoh = 0.15. For comparison, we also 
show the LDOS far away from the impurity site (black line 
in Fig.|^), which resembles the DOS in a pristine bulk. For 
J = 0, since the local hopping is still finite, two in-gap peaks 
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FIG. 5: DFT-based local density of states, a, The local density of states on the IFI, one of its nearest-neighboring Fe sites, that far away 
into the pristine FeSe for a nonmagnetic (upper panel) and spin-polarized IFI (lower panel), b. The d-orbital-resolved partial density of states 
on the nonmagnetic (upper panel) and spin-polarized (lower panel) IFI site. The positive/negative value of each lower panel describes the 
up/down spin component. 


appears near the gap edges. With the increased J, these in-gap 
peaks are shifted toward the Fermi energy (e.g., in the range 
of J G [0,0.8] for ^ = 2.0 as shown in Fig. [^), and cross 
the Fermi energy (e.g., at J ^ 0.8 for ^ = 2.0 as shown in 
Fig.|3^), and move toward the gap edges (e.g., in the range 
of J G [0.8, 2.0] for ^ = 2.0 as shown in Fig. |^). We also 
found a regime of J G [4, 5], for which the in-gap quasipar¬ 
ticle peaks show up again. Since in this large J regime, the 
effect of the topological order becomes non-essential, our re¬ 
sults are comparable with those in the early study on a mag¬ 
netic impurity in an iron-based superconductor |[22l . There¬ 
fore, we do not show here the results for this regime. The in¬ 
gap bound states can also be visualized in the band structure 
in the reduced Brillouin zone of the supercell lattice EB. In 
Fig.[^, all the sub-figures are plotted with varying J for fixed 
values of ^ = 2, Xaoh = 0.15 and A = 0.02. The localiza¬ 
tion nature of these in-gap states are investigated by different 
size of the supercell. As shown in Fig.|^, the dispersiveness 
of the in-gap states is suppressed with increased size of the 
supercell, suggesting that these in-gap states are truly bound 
states. To further investigate the consequence of the topolog¬ 
ical order, we show in Fig.[^ the LDOS on the impurity site 
in the absence (Xaoh = 0) and presence (Xaoh = 0.15) 
of the topological order. The zero-energy peak in the LDOS 
in the absence of the topological order is much narrower (ac¬ 
companied with higher intensity) than that in the presence of 
the topological order, suggesting a much shorter decay length 
of quasiparticle bound state in the former. This contrast is 
further confirmed by the spatial dependence of the LDOS at 
zero-energy, as shown in Fig.|^. Specifically, the zero-energy 
bound state has a longer decay length along the bond direction 
in a two-dimensional lattice. This difference can be under¬ 
stood from the band structure of the pristine system in an un¬ 
folded Brillouin zone, as shown in Fig.|^ where one can find 
that the topological order reduces significantly the gap for the 
wave vector near the X" point along the L — X" — M" path. 


III. RELEVANCE TO RECENT STM EXPERIMENT 

We now discuss the relevance of our results to the recent 
STM experiment on Fe(Te,Se) compound for an adsorbed Fe 
impurity lO. On the one hand, the Fe(Te,Se) compound as 
mentioned in Ref.[3]has the number of (i-electrons correspond¬ 
ing to the half-filling in the two-orbital model, which is differ¬ 
ent from the highly electron-doped KFe 2 Se 2 . Therefore, the 
chemical potential should be set to ensure the half-filling in 
our calculations. In addition, since the distance from the ad¬ 
sorbed Fe site to its nearest neighbors are shorter than the Fe- 
Fe bonding distance of the parent compound, we anticipate 
that the local hopping parameter (0 between the adsorbed 
Fe atom and its neighboring Fe atoms is larger than those 
for the parent compound. On the other hand, since the su¬ 
perconductivity for iron-based compounds arises from the Fe 
3(i-electrons on an essentially square lattice, the normal-state 
part of the model Hamiltonian originally designed for doped 
BaFe 2 As 2 O, to match the low-energy band structure from 
the LDA calculations 12^ . should also be applicable to the 
Fe(Te,Se) system. 

For the present purpose, we carried out the DFT calcula¬ 
tions for an Fe atom adsorbed in the center of Fe plaquette of 
a large real-space structure of Fe(Te,Se) compound. See the 
Methods Section for the calculation details. We have found 
that a spin-polarized Fe atom adsorbed on the surface of an 
otherwise paramagnetic FeSe compound has a lower energy 
than a non-spin-polarized adsorbed Fe atom. Figureshows 
the calculation of LDOS for both non-spin-polarized (upper 
panel) and spin-polarized (lower panel) adsorbed Fe atom. 
The broadened peak of the LDOS on the adsorbed Fe atom 
itself and the smearing of LDOS structure on its neighbor¬ 
ing Fe sites in the FeSe compound (as compared with that 
on the Fe site far way from the adsorbed Fe atom) suggests 
a strong coupling between the adsorbed Fe and its neighbor¬ 
ing Fe sites. It supports the treatment of the local hopping 
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parameters in our model calculations. For the case of the 
spin-polarized case, the magnetization is mostly localized on 
the adsorbed Fe itself with a proximately induced moment on 
its nearest neighboring Fe sites negligibly small. Finally, as 
shown in Fig. Sj, each d-orbital is spin polarized and con¬ 
tributes to the total magnetic moment of 2.89 fiB. However, 
among the total magnetic moment, only about 14.48% (about 
0.42 jib) arises from each of the dxz- and -orbitals, which 
are relevant to the corresponding orbitals active in the pristine 
compound. These finding supports the use of a local magnetic 
impurity with a small spin polarization in our model calcula¬ 
tions. 


IV. DISCUSSION 

Our model study has shown the combined effect of the Z 2 
topological mirror order and the fully gapped s±-wave SC. 
This Z 2 argument extends the mirror symmetry in the mirror- 
Chern representation IT^ for the topological superconductiv¬ 
ity of a 2D lattice. The coexistence of this topological and 
the superconducting order relies on the intrinsic band struc¬ 
ture of the dxz.yz orbitals, which gives a topological metal¬ 
lic state. There the emergence of gapless edge states in the 
coexisting phase is a direct evidence for the topological or¬ 
der. In addition, the introduction of a spin polarized impurity 
can be regarded as a zero-dimension boundary to the 2D lat¬ 
tice, and it induces a similar quasiparticle behavior as around 
a vortex core center in the topological superconductor 1241 . In 
addition, our model calculations of the local electronic struc¬ 
ture have revealed the topological feature of a superconductor 
through the LDOS peak intensity and the size of the impurity- 
induced in-gap bound states. Finally, our DFT calculations of 
an adsorbed Fe atom on the surface of the FeSe compound has 
indeed shown the adsorbed impurity is spin polarized, making 
the findings in the present work relevant to the STM experi¬ 
ment on the Fe(Te,Se) compound. 

Note added. As we nearly completed our research, we 
noticed a recent study by Zhang 1^ that described the in¬ 
gap bound state around an adsorbed Fe atom on the surface 
of iron-based superconductors, within an Anderson impurity 
model with non-spin polarized d-electrons on the impurity 
site. Due to the mismatch of energy levels between the impu¬ 
rity and the Fe-sites in the background, the qualitative agree¬ 
ment with the STM experiment lO was obtained by adjusting 
the chemical potential in the pristine system in the strict con¬ 
dition of an unrealistically small hybridization between the 
adsorbed Fe atom and its neighbors. The shifted chemical po¬ 
tential leads to the system in the highly electron doped regime, 
which is inconsistent with the fact that Fe(Se,Te) is in the 
half-filled regime in the two-orbital model. In addition, it can 
be easily understood that with such a small hybridization, the 
zero-energy peak at the adsorbed Fe site has the character of 
the energy levels of an isolated atom so that the peak intensity 
is too high. 


V. METHODS 

Model Setup. Here we write down the details of = 

T[ti_6, m] + V[\aoh] + 

'T = 'Z 4as (*/J “ M 

I J,a/3,s 

'^= vff c\^^Cjas, ( 3 ) 

{IJ),a,s 

'P= 

IJ,a,s 

were I, J are site index, a, /3 G [dxz^ dyz] are orbital index, 
s G [t? i] is the spin index and /i is the chemical potential 
to adjust the Fermi surface for hall-filling (n = 2). Here, 
G ti_6, are the hopping integrals GJITtI. In our previous 
paper m ^AOH could be obtained self-consistently from the 
NN inter-orbital Coulomb interaction at the mean-field level. 
However, in this work, we follow the same spirit as ref.|26]and 
treat Xaoh and superconducting pair potential A as input pa¬ 
rameters for our analysis. The tensor elements ufj G [0, ±1] 
describe the direction of the NN inter-orbital currents with 

The impurity model, is writ- 

ten, 

■j-IFI ^ ^ tsaicl^a,s <^S,a,s + H.c), 

6,a,s 

~ (^/ 2 ) (4) 

a,ss' 

= (t//2) ^ '^o,Ck:,4.)* 

a 

The index, o, denotes the impurity site and S indicates the NN 
sites from o to its NN sites, tsa is the local hopping integrals 
from the IFI site to its NN sites. Based on the dxz.yz orbital 
orientation, as shown in Fig. [^, there are only two different 
values in our impurity model tsa C here we take to- = 

at^ = ^ with a = —3. In our calculation, different sign and 
magnitude of a do not change our results qualitatively. 

The LDOS can be calculated according to, 

n,k 

k k 

where the and are the eigenfunctions of the BdG 
matrix of the entire Hamiltonian. We use the broadening fa- 
cor F = 0.001 in 5{x) = T + F^). The supercell tech¬ 
nique ED is used for M-repeated cell blocks. 

DFT Calculations. We use the Vienna Ab-initio Simula¬ 
tion Package (VASP) 1^ to carry out the LDA calculations 
for an interstitial Fe impurity (IFI) atom sitting in the center 
of a large real-space structure of Fe(Te,Se) compound. The 
projector augmented planewave method E3 and the Perdew- 
Burke-Emzerhof 1^ exchange-correlation functional were 
adopted. We used a 4 x 4 x 2 Monkhorst-Pack k-point mesh. 
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FIG. 6: Setup of a two-layered with an adsorbed Fe impurity, a, The side/top view of the real-space geometry. The labels 1/2/3 indicate 
the IFI, NN and far-away sites, b, Schematic picture showing the DFT calculated magnetic distribution (black arrow). The actual magnitude 
of each of them is: 2.89/j.b (IFI site), 0.23/j.b (NN site) and 0.019/is (NNN site). In each panel the red/yellow balls represent Fe/Se atoms 
while the blue one stands for the IFI. 


and a 500 eV cutoff. The 2 :-coordinate of Fe impurity atoms 
was fully relaxed up to 10“^ eV/A with other atoms fixed. 

We have considered the setup as shown in Fig.[^, in which 
an IFI atom sitting on the top of a two-layered FeSe. A va¬ 
cancy space in the z-direction was made to mimic the sur¬ 
face effect of the single IFI atom. We have calculated(i) non- 
magnetic(NM) state and (ii) ferro-magnetic(FM) state on the 
IFI site by turning on and off spin polarization initially on 
IFI site. Each of the two calculations has reached the charge 
convergence. Fig.[^ shows the local distribution of magnetic 
moment. The calculated total energy shows that the case (ii) is 
more stable than the case (i), which confirms a spin polarized 
IFI site. 
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